REPORT 1344 


A SIMPLIFIED METHOD FOR APPROXIMATING THE TRANSIENT MOTION IN ANGLES OF 
ATTACK AND SIDESLIP DURING A CONSTANT ROLLING MANEUVER 1 

By Leonabd Stebotield 


SUMMARY 

The transient motion in angles of attack and sideslip during a 
constant rolling maneuver has been analyzed. Simplified ex- 
pressions are presented for the determination of the pertinent 
modes of motion as well as the modal coefficient corresponding 
to each mode. 

Calculations made with and without the derivatives C T p ( side 
force due to sideslip) and C La (lift-curve slope) indicate that 
although these derivatives increase the total damping of the system 
they do not markedly affect the transient motion. 

INTRODUCTION 

Recent flight tests of airplanes designed with their mass 
concentrated primarily in the fuselage have indicated that 
during a rolling maneuver the airplane experiences large 
uncontrollable motions in angles of attack and sideslip. 
A fundamental analysis of this pitch-yaw divergence prob- 
lem, in which the rolling velocity of the airplane is assumed 
constant, is presented in reference 1. This analysis, con- 
cerned primarily with the stability of the system, makes 
possible the calculation of the divergence boundaries and the 
prediction of the range of rolling velocities for which the 
airplane motions will diverge. Analog studies made at the 
National Advisory Committee for Aeronautics for several 
research airplanes have indicated that reference 1 is a helpful 
guide in a roll-coupling investigation, but detailed-motion 
studies based on five degrees of freedom and with pilot inputs 
taken into account are essential to the analysis. Also, 
calculated time histories of the airplane motion indicate that 
largo angles of attack and sideslip which are objectionable to 
a pilot and which may induce severe loads on the airplane 
may be encountered for constant rolling velocities outside 
the critical range. 

The purpose of this paper is to extend the analysis of 
reference 1 by analyzing the transient motion in angles of 
attack and sideslip. Simplified expressions are presented 
which permit an accurate and rapid estimate of the maximum 
angles of attack and sideslip for an airplane rolling at a 
constant velocity which should apply to the case of an air- 
plane entering a rolling maneuver up to the point of recovery 
from the maneuver. The airplane is assumed to be initially 
disturbed by an input in the side-force equation equal to the 
product of the rolling velocity and trim angle of attack. A 
comparison is made of the time histories obtained by using 


the simplified and exact expressions (obtained from the four 
degrees of freedom assuming constant rolling velocity) for 
a currently designed swept-wing fighter flying at a Mach 
number of 0.7 at an altitude of 32,000 feet. 
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wing span, ft 

wing mean aerodynamic chord, ft 


lift coefficient, 


Lift 

ZS 


... , . Rolling moment 

rolling-moment coefficient, 


pitching-moment coefficient, 


Pitching moment 
qSc 


tc • Yawing moment 
yawing-moment coefficient, ^gg 


. , - a. . , Side force 

side-force coefficient, ~ 

gS 

moment of inertia of airplane about principal 
Y-axis, slug-ft 2 

moment of inertia of airplane about principal 
Y-axis, slug-ft 2 

moment of inertia of airplane about principal 
Z-axis, slug-ft 2 
rolling moment, ft-lb 
airplane mass, slugs 
period, sec 

rolling velocity, radians/sec 
steady rolling velocity, radians/sec 
pitching velocity, radians/sec; dynamic pressure, 
lb/sq ft 

yawing velocity, radians/sec 
wing area, sq ft 
time, sec 

time to damp to one-half amplitude, sec 
time to double amplitude, sec 
airplane velocity, ft/sec 

initial angle of attack of airplane principal axis, 
radians 

incremental change in angle of attack, radians 


a=ct (,-j-Aa 


£ angle of sideslip, radians 

X root of characteristic equation 

adzia complex roots of characteristic equation 

AjjXj real roots of characteristic equation 


i Supersedes recently declassified NACA Research Memorandum L66F04 by Leonard Stemfleld^IGW 
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A a= (Aa)^ 1 are used in the equations written in determinant 
form, X must be a root of the stability equation 
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EQUATIONS OF MOTION 


In this analysis, on the assumption of constant forward 
velocity and constant rolling velocity, the following linearized 
equations of motion of an airplane referred to the principal 
axes were used: 

Pitching: 



i ( Iz—Ix 

\ M.q 

. M a 


2 \ It , 

fj? 
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Tawing: 





- (Ix-E 

Nr 

.Nf, 


\ Iz 


r+ ~r 

-Lz 

Side force: 





P+r—po 

Aa =(^f-P+P° a * 

Normal force: 





i-q+&P=-^f 

■ r Aa 


(la) 

(lb) 

(lc) 

(l d) 


When the substitutions q=q<A i , r=r 0 e M , j3=)3 0 e x ‘, and 


L g Yp ( lx Iy\ ( Iz Ix \ 

mVmV\ I z A It ) 


V 2 


In reference 1, G L<t and 0 Y p were assumed to be equal to 
zero; this assumption is representative of the condition in 
which the center of gravity is moving in a straight lino. In 
this analysis the case where C La =C Y p~d (referred to as case 
(a)) is considered, as well as the case where the values of 
C La and 0 Y p are finite (referred to as case (b)). 

ANALYSIS AND DISCUSSION 

DIVERGENCE BOUNDARIES 

The conditions necessary for the system described by 
equations (1) to be stable are that the coefficients of the 
characteristic equation (eq. (2)) must be positive and that 
South's discriminant, BCD—AIP—IPE, must be greater 
than zero. From an examination of the characteristic equa- 
tion of this system, it can be shown that, for aircraft designed 
with positive static stability and for conditions where each 
of the damping derivatives Onj, O m g, 0 La , and 0 Y& adds damp- 
ing to the system, the oscillatory modes will be stable for all 
values of the steady rolling velocity. Thus, the only type of 
instability that could be encountered is an aperiodic mode 
which occurs when the constant term of the characteristic 
equation, the E coefficient, is negative. As shown in ref- 
erence!, divergence boundaries, obtained by setting E= 0, 
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can bo plotted as a function of the squares of the natural 
frequencies in pitch and yaw nondimensionalized to the 
square of the rolling velocity, cor and co/, respectively. 
These divergence boundaries are shown in figure 1 for cases 
(a) and (b) for an airplane whose mass and aerodynamic 
characteristics are presented in table I. The nondimensional 
pitch and yaw frequencies of the airplane described in 
table X fall along the straight line, with the slope equal 
^ J[ ~c, 

2 shown in figure 1. Each point on this line 


to 




<-Z _ 

Jr b’ 


corresponds to a particular value of p a and, as p< s increases, 
the point moves along the line toward the origin. For case 
(a), the airplane line intersects the divergence boundaries 
at two points and thus defines the critical range of roll- 
ing velocities where a divergent mode exists, namely, 
1.86<|poK2.33. For case (b) the airplane line lies in the 


TABLE I 

MASS AND AERODYNAMIC CHARACTERISTICS OF 
AIRPLANE 


m, slugs 745 

lx, slug-ft 3 10, 976 

1 y, slug-ft 3 _ 57,100 

lz, slug-ft 3 64,975 

C„ T , per radian - —0.095 

C„ a , per radian -. —3.6 

C„ a , per radian _ —0. 36 

C „0 , per radian — - 0. 067 

Cr 0 , per radian- —0. 28 

Cl„, per radian- 3. 85 

Ci v , per radian,. - - —0. 255 

F, ft/sec.. 691 

5, sq ft 1 — 377 

c, ft 11.3 

6, ft 3a 6 

q, lb/sq ft 197 



Figure 1. — Boundaries in the wp, plane which define Jregions of 
aperiodic divergence for example airplane. 


stable region and the motion is therefore stable for all 
values of p 0 . 

CHARACTERISTIC MODES 

The roots of the characteristic stability equation cor- 
responding to the absolute values of p 0 shown in figure 1 are 
presented in tables II (a) and II (b) for cases (a) and (b), 


TABLE n 

ROOTS OF THE CHARACTERISTIC EQUATION FOR VALUES OF p„ SHOWN IN FIGURE 1 

(a) Roots for case (a) 


Vo 

Exact 

Approximate 

0 

-0. 210 ±2. 29* 

-0. 0526 ±1. 54* 



-L 0 

—0. 156 ±2. 90* 

-0. 107 + 0. 922* 

-0. 152+3. 06* 

-0. 111 + a 874* 

-1. 5 

-0. 143 ±3. 34* 

-0. 12 ±0. 464* 

-0. 144+3. 39* 

-0. 119 + 0. 457* 

— 1. 86 

-0. 137 ±3. 66* 

-0. 251 0 

-0. 133 + a 67* 

-0. 25 0 

-2. 0 

-0. 135 ±3. 79* 

-0. 355+0. 0996 

-0. 136 + 3. 79 

-0. 354+0. 0997 

-2. 33 

-0. 131 ±4. 09* 

-0. 256 0 

-0. 132+4. 10* 

-0. 252 0 

-2.6 

-0. 129 ±4. 24* 

-0. 134+0. 267* 

-0. 130+4. 26* 

-a 133 + 0. 294*' 

-3. 0 

-0. 124 ±4. 70* 

-0. 139 + 0. 768* 

-0. 125 + 4. 77* 

-0. 138 + a 756* 


(b) Roots for case (b) 


Po 

- 

Exact 

Approximate 

0 

-L 0 

' -0.488 +2.30* 
-0. 362+2. 89* 

-0. 0729 + 1. 54* 
-0. 199 ±0. 942* 

- 0. 352 ± 3. 10* - 0. 209+ 0. 876* 

-1. 5 

-0. 337+3. 33* 

-0.224 ± 0.483* 

- 0. 340 ± 3. 42* - 0. 221 ± a 469* 

-L 86 

-0. 327+3. 66* 

-0. 322 -0. 145 

-0.332 + 3.71* -0.311-0.146 

-2. 0 

-0. 324+3. 79* 

— 0. 453 —0. 020 

- 0. 329 ± 3. 83* - 0. 444- 0. 020 

-2. 33 

-0. 318+4, 08* 

-0. 374 -0. Ill 

-0. 322++ 13* -0. 363 - 0. 113 

-2.5 

-0. 316 + 4, 24* 

-0.245 ±0.258* 

-0. 320 + + 29* -0.241 + 0.255* 

-3. 0 

-0. 311 + 4. 70* 

-0. 250 ±0. 760* 

— 0I312++ 80* -0.248 + 0.743* 
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respectively. A method of approximating the roots of the 
characteristic equation is presented in the appendix. The 
roots for p s =0 which correspond to the condition where no 
coupling exists between the lateral and longitudinal modes 
are also given in the table. If the roots are complex, 

\=a±ioi, P=—> and Ty— — whereas, for a real root 


Xi> 


T }i = 


—0.693 

Xi 


or T 2 = 


0.693 

Xt 


The roots for p o =0, for both eases (a) and (b), indicate that 
the motion consists of two independent oscillatory modes, a 
longitudinal and lateral oscillation with periods of 2.74 
seconds and 4.08 seconds, respectively. The damping of the 
oscillation for case (b) is greater than that for case (a) because 
the damping derivatives Q r and C La were assumed equal to 
zero in case (a). The values of Ty in seconds are given in 
the following table: 




Cane 




Longitudinal 

Lateral 


oscillation 

oscillation- 


3. 3 

13. 2 

■■ 

L 4 

9. 5 


For values of [pol other than zero, the lateral variables /3 
and r are coupled with the longitudinal variables a and q 
and the system now has the characteristic modes shown in 
table II for each of the variables j3, r, a, and q. It should be 
noted that as |p 0 | increases, the frequency of the high fre- 
quency mode continues to increase and the damping de- 
creases. However, the frequency of the other oscillation 
decreases until |p 0 | reaches approximately 1.86 when the 
oscillation breaks down into two aperiodic modes. In case 
(a) at this value of |p 0 |, a zero root is obtained and the un- 
stable region begins. The system remains unstable, with 
one aperiodic divergent mode, until \p 0 \ becomes equal to 
2.33. In case (b), the system is stable although two aperiodic 
modes occur for 1.86<|p 0 K2.33. For [p 0 |>2.33, the oscilla- 
tion reappears in cases (a) and (b) and the frequency and 
damping increase as |p 0 | increases. 

From the contour lines of constant oscillation frequencies 
presented in reference 1, it appears that, in general, the 
frequency of one of the oscillatory modes will always decrease 
and probably break down into two aperiodic modes as [p 0 | 
increases. 


MOTIONS IN (3(f) AND Ao<f) 

Time histories were calculated by the method of Laplace 
transform (ref. 2) of the motions in and Aa for the several 
values of p Q shown in figure 1. The disturbance acting on 
the airplane is represented by the term p 0 a « in the side-force 
equation while the airplane is rolling in a negative direction. 

The motions in and a — are presented in figures 2 and 
ao a 0 

3 for cases (a) and (b), respectively. A comparison of figures 
2 and 3 shows that the inclusion of C\ a and modifies 




(b) p 0 “= — 1-5. 

Figure 2. — Time histories of 0/ao and Aa/ao for ease (a) . 

the motions slightly. Additional time histories were made 
and are presented in figure 4 for the same airplano but 
with the following assumed values of the derivatives: 
Om a =— 0.09, 0^=0.114, and Cy e =C La =0. Forp 0 = — 1.0, 
the airplane would be located in figure 1 at the point 
o> a 2 = 1.32 and «* 2 = 4.77. 

Although the modes of motion involved in the transient 
behavior of the airplane will be stable for points along 
the airplane line located in the stable region, that is, 
2.33<C|j l ol'0-86 in case (a) and for all values of |p 0 | in case (b), 
the magnitudes of Aa and /S developed during the tran- 
sient motion may be large enough to cause severe loads on 
the airplane. For example, in the particular flight condition 
for which the motions were calculated, where a 0 =5°, the 
airplane develops a sideslip angle of about —8° at po= — 1.5 
and Aa is about —20° at po= — 3.0. These values are 
reached in about 3 seconds. 

In all the motions the high-frequency mode is clearly evi- 
dent, although the amplitude of this mode is relatively small. 
Hence, the high-frequency mode modifies the shape of the 
motion only slightly. Thus, for p 0 = — 1.0, —1.5, —2.5, 
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(e) po= — 2.33. 

Figubb 2. — Continued. 

and —3.0, the low-frequency mode is of primary impor- 
tance in the motions, whereas for p 0 = — 1.86, —2.0, and 
— 2.33, the aperiodic modes are of primary importance. 
In addition, the analytical expressions for the motions 
would also contain a constant term which represents the 
stdady-state value of (3 and A a if the system is stable. A 
good approximation to the time history could therefore be 
expressed in general form for p 0 = — 1-0, —1.5, —2.5, and 
—3.0 as 

M or ^H =Ko+K e ai sin ( wf+e ) (3) 

Oo <*0 

For case (b) when p 0 = — 1.86, —2.0, and —2.33 and for 
case (a) when p 0 =— 2.0, the following equation is obtained: 

M or ^.^Kc+K^+K,^ 1 

Oo «0 


(d) p 0 =— 2.0. 
Figube 2. — Continued. 


( 4 ) 
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Figure 2. — Concluded. 


I 


Figure 3. — Time histories of /J/ao and Aajao for ease (b). 









3 4 

/, sec 


(e) pa—— 2.33. 
Figure 3. — Continued. 


3 4 5 6 7 

/, sec 

(d) p 0 ~-2.0. 

Figube 3. — Continued. 

and for case (a) when p 0 = — 1.86 and —2.33, 

^or^^+Jf^-1) (5) 

0(q CtQ 

Since the high-frequency mode is omitted from expressions 

(3), (4), and (5), the values of and — a ^ will not satisfy 

<*0 ®o 

the initial condition of being zero at f=0. 

By differentiating expressions (3), (4), and (5) and setting 
them equal to zero, the maximum values of /3 and Aa and the time at which they occur could be easily determined. In these 
calculations, /3 and Aa occurred at about J4 and ji of the period, respectively. 

APPROXIMATION OF THE AMPLITUDE COEFFICIENTS AND PHASE ANGLE 

The amplitude coefficients Ko, K lt Ki, and K 6 and the phase angle t in expressions (3), (4), and (5) can be calcu- 

lated directly by using the Heaviside expansion formula given on page 45 of reference 2. A detailed examination of the 
exact expicssions for the amplitude coefficients and phase angle indicated that many of the terms appearing in the expressions 
had a negligible effect on the resultant magnitudes. Thus the following simplified expressions were derived which result in 
good agreement with the exact values. 

In the expressions for 
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Kj,K 3 , and K, 




2 C\+D 


In the expressions for ■ 


*-=g {f f+Mf t, «(¥)(¥)]} 


k x =JgFFc? 

A 2) i£ 3 , and its— - X (2C r X+D) 


*•=-£ [(*£*) » --f (^>f f ] 


c^-cA-. 


(a\ *■-»’(%+%) 

\w/ to (C—2u?) 


—tan -1 — The (7, Z>, and E appearing in the preceding expressions 

6 Ci are the coefficients of the characteristic equation. A com- 

w h ere parison of the motions obtained by using the exact and 

approximate coefficients is shown in figures 2 to 4 to be very 

c _ e? ^ _|7 kzM) W-p* ( Iz ~ —} ( g° od - 

2 CO 5 (C — 2 to 2 ) |_ V Iy J Iz V I? J \ lz / J EFFECT OF ASSUMING po CONSTANT AT t — 0 

^ In general, the rolling velocity builds up exponentially in 

\ response to an aileron input instead of being constant at 

\ — £_q ^ aggujned this analysis. In order to determine the 

_i effect of the assumption on the motions in j8 and Aa, an 

' X analog study was conducted with the assumption that the 

N '\^ " rolling velocity reaches its steady-state value exponentially. 

_ 2 The resultant motions in /3 and Aa were compared with the 





Approximate expression 

Exact expression 



uu mmm 

min 


0 I 


2 3 4 

/, sec 

(f) po= — 2.6. 
Figure 3. — Continued. 


(g) po=— 3.0. 
Figure 3. — Concluded. 
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Figure 4. — Time histories of p/a 0 and A«/a 0 for C„ a = —0.09, 
C* f — 0.114, po= — 1.0, and Cy fi =C La = 0. 


motions in /S 'and A a obtained on the assumption that p 0 is 
constant at t=Q. Equations (1) wore used in the analog 
study and, in addition, a simplified roll equation 

was introduced. This equation states that 

p reaches its steady-state value, p a , according to the expression 
p=Po(.e u —l). In general, the resultant motions obtained 
in the analog study were very similar to the motions shown 
in figures 2 and 3. The only significant difference noted 

was that d i fferent peak values of — and— were encountered. 

«0 Oo 

Table III presents a comparison between the maximum 

values of — • and — obtained from the analog study and 
«o «o 

from figures 2 and 3 for p Q — —1.0, —1.5, and —3.0. It is 
seen from the values presented in table III that higher 
maximum values occur for p 0 - —1.0 and —1.5 for the con- 


dition of constant p 0 at f=0 whereas the opposite is true for 
Po = — 3.0. These results can be explained by the fact that, 
during the first several seconds of the transient motion, the 
average value of p 0 in the analog study is smaller than under 
the assumption of constant p a . Hence in the analog study 
for po = —1.0 and —1.5 the airplane may be considered to 
be located farther away from the stability boundary, but 
for p 0 = —3.0 the airplane is located closer to the stability 
boundary. 


CALCULATIONS OP MOTIONS FOS AN AIRPLANE PERFORMING A 
360° ROLL 


The analysis presented thus far is applicable to the con- 
dition of the airplane performing a continuous rolling motion. 
Of particular interest are the motions of the airplane when 
the pilot performs a 360° roll. These motions can be 
approximated by first determining the values of fi, r, a , and 
q at the time the airplane has rolled through 360°, that is, 
2 T 

at t=— At this value of t, the constant rolling velocity 

is returned to p 0 =0 and the motion originally described 
by four degrees of freedom is now separable into its lateral 
and longitudinal values, the former represented by the 
yawing and sideslipping equations and the latter repre- 
sented by the pitching-moment and normal-force equations. 


Thus, the values of r, a, and q at £=— are the initial 

p 0 

conditions required to calculate the motions subsequent to 

2 tt 

t = — . The values of and a can be determined from 
Po 

equation (3), (4), or (5). From equations (lc) and (Id), 
the following equations are obtained: 


and 


7, 


r=^y-/3+2> o a°+Po Aa— 


2=“|r Aa +?W3+« 


( 6 ) 

( 7 ) 


2v 

A good approximation of r and q at f=— is obtained by 

neglecting the terms /3 and a in equations (6) and (7). The 

2tt 

following table compares the values of r and q at t= — for 

Po 

p 0 =—1.5, —1.7, and —3.0 when and a are included in 


TABLE III 

COMPARISON OF MAXIMUM VALUES OF - AND — 

; CCq Oto 


Po 

C'v Cy fi 

( £ ) 

VRl /max 

/ Aa 
\«o 

^ nai 

Po constant at (=0 

•e 

II 

S 

% 

\ 

t— * 

Po constant at t= 0 

P=P o (e*‘— 1) 

-L 0 

0 

-0. 64 

-0. 48 

0. 6 

0.4 

-L 0 

Finite 

-. 6 

-. 48 

. 4 

. 3 

-L 5 

0 

-L 64 

-1. 27 

L 0 

. 9 \ 

-1. 6 

Finite 

-L 48 

-L 12 

. 68 

. 6 

-3.0 

0 

— L 75 

-L 91 

-3. 92 

—4. 15 

-3.0 

Finite 

-L 90 

-2. 1 

-a 25 

- -3.3 
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equations (6) and (7) and then deleted from equations (6) 
and (7): 


Po 

0 and h included 

jS and a deleted 

r 

2 

r 

2 

-1.6 

-a 20 

0. 20 

-0. 19 

0. 19 

-1. 7 

-. 12 

. 33 

-. 14 

.32 

-3. 0 

. 26 

. 30 

. 29 

.40 


EFFECT OF GRAVITY TERMS ON THE P AND Ar MOTIONS 


The effect of the gravity components which would appear 
in the side force (eq. (lc)) and normal force (eq. (Id)) equa- 
tions have been neglected in this analysis. In order to de- 
termine the effect of the gravity terms on the and Aa 
motions, an analog study was conducted using equations 
(lc) and (Iff) with the gravity terms included and then 
deleted from the equations. The airplane was assumed to 
be performing a 360° roll for rolling velocities of %><,=— 1.5, 

— 1.7, and —3.0. A comparison of the maximum values of 

— and — obtained with and without the gravity terms is 
«o ao 

presented in table IV. The two values given for that 


is, and correspond to the maximum values of — 
ao ao <*o 

obtained during the initial part of the transient motion 

(pa= constant) and during the recovery part of the motion 

(p 0 =0), respectively. The comparison shows that with the . 

gravity terms included the maximum values of /S are about 

2° greater and the maximum values of Aa increase by about 

1° for an initial ao of 5°. 


TABLE IV 

COMPARISON OF MAXIMUM VALUES OF - AND — 

ao ao 


Vo 

Gravity terms included 

Gravity terms deleted 

(“) 
\&o/ max 

(-) 
\«o /max 

(~) 

\ CCO /max 

( £ ) 

\ao/ma* 

(— ) 

\ «0 /max 

/AaA 

\ O0 ) max 

— L 6 
-h T 
-3. 0 

— L 8 
-2. 5 
-3.0 

—0. 50 
-. 60 
-2. 4 

0. 90 
L 3 

1. 6 

-1.4 
-2. 1 
-2.7 

-0. 30 
-.40 
-2. 2 

H 


CONCLUDING REMARKS 

From the analysis presented in this paper it appears that 
the transient motion in angles of attack and sideslip during a 
constant rolling maneuver consists chiefly of either an oscilla- 
tory mode or two aperiodic modes. Approximate expressions 
are derived for the determination of these modes as well as 
the modal coefficient corresponding to each mode. Inclusion 
of the derivatives 0 Yp (side force due to sideslip) and Ol b 
(lift-curve slope) increases the total damping of the system 
but does not markedly affect the transient motions. 

The sole input considered in this paper is the term p f fit a 
(the product of rolling velocity and trim angle of attack) in 
the side-force equation. 
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APPENDIX 

METHOD OP APPROXIMATING THE ROOTS OP THE CHARACTERISTIC EQUATION 


Tho fourth-degree characteristic equation 
AX 4 +5X 8 + O\ 3 +D\+E=0 
may be factored as follows: 

where 

A= 1 

B=ai+a^ 

0=biAbi~{~(Li(h 

T) =z (L 3 b\~\“Q / ib} 

E=bib 3 

An examination of tho quadratic equations from which the 
exact roots presented in tables II (a) and II (b) were obtained 
indicated that b{^b 3 and that a L and a% are of the same order 
of magnitude and much smaller than b\. Thus, one may 
write 

A= 1 

jB=a 1 +a 3 


C=b, 

D=aJ>i 

E=b 1 b 1 

The solution of these equations yields 

BO-D 
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b x =C 

D 
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The roots of X s +OiX+6i=0 approximate the high-fre- 
quency mode whereas the roots of X 2 +d2X+5 2 =0 approxi- 
mate the remaining two roots of primary importance in the 
motion calculations presented herein. A comparison between 
the exact roots and the approximate roots presented in tables 
H (a) and II (b) indicates that the agreement is excellent. 
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